Calling Structural Variants with Confidence from Short-Read Data in Wild Bird Populations.

Comprehensive characterization of structural variation in natural populations has only become feasible in the last decade. To investigate the population genomic nature of structural variation, reproducible and high-confidence structural variation callsets are first required. We created a population-scale reference of the genome-wide landscape of structural variation across 33 Nordic house sparrows (Passer domesticus). To produce a consensus callset across all samples using short-read data, we compare heuristic-based quality filtering and visual curation (Samplot/PlotCritic and Samplot-ML) approaches. We demonstrate that curation of structural variants is important for reducing putative false positives and that the time invested in this step outweighs the potential costs of analyzing short-read-discovered structural variation data sets that include many potential false positives. We find that even a lenient manual curation strategy (e.g. applied by a single curator) can reduce the proportion of putative false positives by up to 80%, thus enriching the proportion of high-confidence variants. Crucially, in applying a lenient manual curation strategy with a single curator, nearly all (>99%) variants rejected as putative false positives were also classified as such by a more stringent curation strategy using three additional curators. Furthermore, variants rejected by manual curation failed to reflect the expected population structure from SNPs, whereas variants passing curation did. Combining heuristic-based quality filtering with rapid manual curation of structural variants in short-read data can therefore become a time- and cost-effective first step for functional and population genomic studies requiring high-confidence structural variation callsets.

The genetic basis of dispersal in a vertebrate metapopulation.

Dispersal affects evolutionary processes by changing population size and genetic composition, influencing the viability and persistence of populations. Investigating which mechanisms underlie variation in dispersal phenotypes and whether populations harbour adaptive potential for dispersal is crucial to understanding the eco-evolutionary dynamics of this important trait. Here, we investigate the genetic architecture of dispersal among successfully recruited individuals in an insular metapopulation of house sparrows. We use an extensive long-term individual-based ecological data set and high-density single-nucleotide polymorphism (SNP) genotypes for over 2500 individuals. We conducted a genome-wide association study (GWAS), and found a relationship between dispersal probability and a SNP located near genes known to regulate circadian rhythm, glycogenesis and exercise performance, among other functions. However, this SNP only explained 3.8% of variance, suggesting that dispersal is a polygenic trait. We then used an animal model to estimate heritable genetic variation (σA 2 ), which composes 10% of the total variation in dispersal probability. Finally, we investigated differences in σA 2 across populations occupying ecologically relevant habitat types (farm vs. non-farm) using a genetic groups animal model. We found different adaptive potentials across habitats, with higher mean breeding value, σA 2 , and heritability for the habitat presenting lower dispersal rates, suggesting also different roles of environmental variation. Our results suggest a complex genetic architecture of dispersal and demonstrate that adaptive potential may be environment dependent in key eco-evolutionary traits. The eco-evolutionary implications of such environment dependence and consequent spatial variation are likely to become ever more important with the increased fragmentation and loss of suitable habitats for many natural populations.

Maternal transmission gives way to social transmission during gut microbiota assembly in wild mice.

BACKGROUND: The mammalian gut microbiota influences a wide array of phenotypes which are relevant to fitness, yet knowledge about the transmission routes by which gut microbes colonise hosts in natural populations remains limited. Here, we use an intensively studied wild population of wood mice (Apodemus sylvaticus) to examine how vertical (maternal) and horizontal (social) transmission routes influence gut microbiota composition throughout life. RESULTS: We identify independent signals of maternal transmission (sharing of taxa between a mother and her offspring) and social transmission (sharing of taxa predicted by the social network), whose relative magnitudes shift as hosts age. In early life, gut microbiota composition is predicted by both maternal and social relationships, but by adulthood the impact of maternal transmission becomes undetectable, leaving only a signal of social transmission. By exploring which taxa drive the maternal transmission signal, we identify a candidate maternally-transmitted bacterial family in wood mice, the Muribaculaceae. CONCLUSION: Overall, our findings point to an ontogenetically shifting transmission landscape in wild mice, with a mother's influence on microbiota composition waning as offspring age, while the relative impact of social contacts grows.

Assortative Mating in an Ecological Context: Effects of Mate Choice Errors and Relative Species Abundance on the Frequency and Asymmetry of Hybridization.

AbstractThe frequency and asymmetry of mixed-species mating set the initial stage for the ecological and evolutionary implications of hybridization. How such patterns of mixed-species mating, in turn, are influenced by the combination of mate choice errors and relative species abundance remains largely unknown. We develop a mathematical model that generates predictions for how relative species abundances and mate choice errors affect hybridization patterns. When mate choice errors are small (<5%), the highest frequency of hybridization occurs when one of the hybridizing species is at low abundance, but when mate choice errors are high (>5%), the highest hybridization frequency occurs when species occur in equal proportions. Furthermore, females of the less abundant species are overrepresented in mixed-species matings. We compare our theoretical predictions with empirical data on naturally hybridizing Ficedula flycatchers and find that hybridization is highest when the two species occur in equal abundance, implying rather high mate choice errors. We discuss ecological and evolutionary implications of our findings and encourage future work on hybrid zone dynamics that take demographic aspects, such as relative species abundance, into account.

Infection by a helminth parasite is associated with changes in DNA methylation in the house sparrow.

Parasites can exert strong selective pressures on their hosts and influence the evolution of host immunity. While several studies have examined the genetic basis for parasite resistance, the role of epigenetics in the immune response to parasites is less understood. Yet, epigenetic modifications, such as changes in DNA methylation, may allow species to respond rapidly to parasite prevalence or virulence. To test the role of DNA methylation in relation to parasite infection, we examined genome-wide DNA methylation before and during infection by a parasitic nematode, Syngamus trachea, in a natural population of house sparrows (Passer domesticus) using reduced representation bisulfite sequencing (RRBS). We found that DNA methylation levels were slightly lower in infected house sparrows, and we identified candidate genes relating to the initial immune response, activation of innate and adaptive immunity, and mucus membrane functional integrity that were differentially methylated between infected and control birds. Subsequently, we used methylation-sensitive high-resolution melting (MS-HRM) analyses to verify the relationship between methylation proportion and S. trachea infection status at two candidate genes in a larger sample dataset. We found that methylation level at NR1D1, but not CLDN22, remained related to infection status and that juvenile recruitment probability was positively related to methylation level at NR1D1. This underscores the importance of performing follow-up studies on candidate genes. Our findings demonstrate that plasticity in the immune response to parasites can be epigenetically mediated and highlight the potential for epigenetic studies in natural populations to provide further mechanistic insight into host-parasite interactions.

Functional genomic tools for emerging model species.

Most studies in the field of ecology and evolution aiming to connect genotype to phenotype rarely validate identified loci using functional tools. Recent developments in RNA interference (RNAi) and clustered regularly interspaced palindromic repeats (CRISPR)-Cas genome editing have dramatically increased the feasibility of functional validation. However, these methods come with specific challenges when applied to emerging model organisms, including limited spatial control of gene silencing, low knock-in efficiencies, and low throughput of functional validation. Moreover, many functional studies to date do not recapitulate ecologically relevant variation, and this limits their scope for deeper insights into evolutionary processes. We therefore argue that increased use of gene editing by allelic replacement through homology-directed repair (HDR) would greatly benefit the field of ecology and evolution.

Effects of environment and genotype on dispersal differ across departure, transfer and settlement in a butterfly metapopulation.

Active dispersal is driven by extrinsic and intrinsic factors at the three stages of departure, transfer and settlement. Most empirical studies capture only one stage of this complex process, and knowledge of how much can be generalized from one stage to another remains unknown. Here we use genetic assignment tests to reconstruct dispersal across 5 years and 232 habitat patches of a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We link individual dispersal events to weather, landscape structure, size and quality of habitat patches, and individual genotype to identify the factors that influence the three stages of dispersal and post-settlement survival. We found that nearly all tested factors strongly affected departure probabilities, but that the same factors explained very little variation in realized dispersal distances. Surprisingly, we found no effect of dispersal distance on post-settlement survival. Rather, survival was influenced by weather conditions, quality of the natal habitat patch, and a strong interaction between genotype and occupancy status of the settled habitat patch, with more mobile genotypes having higher survival as colonists rather than as immigrants. Our work highlights the multi-causality of dispersal and that some dispersal costs can only be understood by considering extrinsic and intrinsic factors and their interaction across the entire dispersal process.

Photoperiod controls wing polyphenism in a water strider independently of insulin receptor signalling.

Insect wing polyphenism has evolved as an adaptation to changing environments and a growing body of research suggests that the nutrient-sensing insulin receptor signalling pathway is a hot spot for the evolution of polyphenisms, as it provides a direct link between growth and available nutrients in the environment. However, little is known about the potential role of insulin receptor signalling in polyphenisms which are controlled by seasonal variation in photoperiod. Here, we demonstrate that wing length polyphenism in the water strider Gerris buenoi is determined by photoperiod and nymphal density, but is not directly affected by nutrient availability. Exposure to a long-day photoperiod is highly inducive of the short-winged morph whereas high nymphal densities moderately promote the development of long wings. Using RNA interference we demonstrate that, unlike in several other species where wing polyphenism is controlled by nutrition, there is no detectable role of insulin receptor signalling in wing morph induction. Our results indicate that the multitude of possible cues that trigger wing polyphenism can be mediated through different genetic pathways and that there are multiple genetic origins to wing polyphenism in insects.

Wild epigenetics: insights from epigenetic studies on natural populations.

Epigenetic mechanisms such as DNA methylation, histone modifications and non-coding RNAs are increasingly targeted in studies of natural populations. Here, I review some of the insights gained from this research, examine some of the methods currently in use and discuss some of the challenges that researchers working on natural populations are likely to face when probing epigenetic mechanisms. While studies supporting the involvement of epigenetic mechanisms in generating phenotypic variation in natural populations are amassing, many of these studies are currently correlative in nature. Thus, while empirical data point to widespread contributions of epigenetic mechanisms in generating phenotypic variation, there are still concerns as to whether epigenetic variation is instead ultimately controlled by genetic variation. Disentangling these two sources of variation will be a key to resolving the debate about the importance of epigenetic mechanisms, and studies on natural populations that partition the relative contribution of genetic and epigenetic factors to phenotypic variation can play an important role in this debate.

Rapid evolution of sexual size dimorphism facilitated by Y-linked genetic variance.

Sexual dimorphism is ubiquitous in nature but its evolution is puzzling given that the mostly shared genome constrains independent evolution in the sexes. Sex differences should result from asymmetries between the sexes in selection or genetic variation but studies investigating both simultaneously are lacking. Here, we combine a quantitative genetic analysis of body size variation, partitioned into autosomal and sex chromosome contributions and ten generations of experimental evolution to dissect the evolution of sexual body size dimorphism in seed beetles (Callosobruchus maculatus) subjected to sexually antagonistic or sex-limited selection. Female additive genetic variance (VA) was primarily linked to autosomes, exhibiting a strong intersexual genetic correlation with males ([Formula: see text] = 0.926), while X- and Y-linked genes further contributed to the male VA and X-linked genes contributed to female dominance variance. Consistent with these estimates, sexual body size dimorphism did not evolve in response to female-limited selection but evolved by 30-50% under male-limited and sexually antagonistic selection. Remarkably, Y-linked variance alone could change dimorphism by 30%, despite the C. maculatus Y chromosome being small and heterochromatic. Our results demonstrate how the potential for sexual dimorphism to evolve depends on both its underlying genetic basis and the nature of sex-specific selection.

Dispersal in a house sparrow metapopulation: An integrative case study of genetic assignment calibrated with ecological data and pedigree information.

Dispersal has a crucial role determining ecoevolutionary dynamics through both gene flow and population size regulation. However, to study dispersal and its consequences, one must distinguish immigrants from residents. Dispersers can be identified using telemetry, capture-mark-recapture (CMR) methods, or genetic assignment methods. All of these methods have disadvantages, such as high costs and substantial field efforts needed for telemetry and CMR surveys, and adequate genetic distance required in genetic assignment. In this study, we used genome-wide 200K Single Nucleotide Polymorphism data and two different genetic assignment approaches (GSI_SIM, Bayesian framework; BONE, network-based estimation) to identify the dispersers in a house sparrow (Passer domesticus) metapopulation sampled over 16 years. Our results showed higher assignment accuracy with BONE. Hence, we proceeded to diagnose potential sources of errors in the assignment results from the BONE method due to variation in levels of interpopulation genetic differentiation, intrapopulation genetic variation and sample size. We show that assignment accuracy is high even at low levels of genetic differentiation and that it increases with the proportion of a population that has been sampled. Finally, we highlight that dispersal studies integrating both ecological and genetic data provide robust assessments of the dispersal patterns in natural populations.

Social networks strongly predict the gut microbiota of wild mice.

The mammalian gut teems with microbes, yet how hosts acquire these symbionts remains poorly understood. Research in primates suggests that microbes can be picked up via social contact, but the role of social interactions in non-group-living species remains underexplored. Here, we use a passive tracking system to collect high resolution spatiotemporal activity data from wild mice (Apodemus sylvaticus). Social network analysis revealed social association strength to be the strongest predictor of microbiota similarity among individuals, controlling for factors including spatial proximity and kinship, which had far smaller or nonsignificant effects. This social effect was limited to interactions involving males (male-male and male-female), implicating sex-dependent behaviours as driving processes. Social network position also predicted microbiota richness, with well-connected individuals having the most diverse microbiotas. Overall, these findings suggest social contact provides a key transmission pathway for gut symbionts even in relatively asocial mammals, that strongly shapes the adult gut microbiota. This work underlines the potential for individuals to pick up beneficial symbionts as well as pathogens from social interactions.

Temporal changes in DNA methylation and RNA expression in a small song bird: within- and between-tissue comparisons.

BACKGROUND: DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled. RESULTS: We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body. CONCLUSION: Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case. However, associations between temporal changes in DNA methylation with changes in gene expression are mostly tissue- and genomic location-dependent. The observation that temporal changes in DNA methylation within RBCs can relate to changes in gene expression in less accessible tissues is important for a better understanding of how environmental conditions shape traits that temporally change in expression in wild populations.

Rapid changes in DNA methylation associated with the initiation of reproduction in a small songbird.

Species with a circannual life cycle need to match the timing of their life history events to the environment to maximize fitness. However, our understanding of how circannual traits such as timing of reproduction are regulated on a molecular level remains limited. Recent studies have implicated that epigenetic mechanisms can be an important part in the processes that regulate circannual traits. Here, we explore the role of DNA methylation in mediating reproductive timing in a seasonally breeding bird species, the great tit (Parus major), using genome-wide DNA methylation data from individual females that were blood sampled repeatedly throughout the breeding season. We demonstrate rapid and directional changes in DNA methylation within the promoter region of several genes, including a key transcription factor (NR5A1) known from earlier studies to be involved in the initiation of timing of reproduction. Interestingly, the observed changes in DNA methylation at NR5A1 identified here are in line with earlier gene expression studies of reproduction in chicken, indicating that the observed shifts in DNA methylation at this gene can have a regulatory role. Our findings provide an important step towards elucidating the genomic mechanism that mediates seasonal timing of a key life history traits and provide support for the idea that epigenetic mechanisms may play an important role in circannual traits.

On the Use of Blood Samples for Measuring DNA Methylation in Ecological Epigenetic Studies.

There is increasing interest in understanding the potential for epigenetic factors to contribute to phenotypic diversity in evolutionary biology. One well studied epigenetic mechanism is DNA methylation, the addition of a methyl group to cytosines, which have the potential to alter gene expression depending on the genomic region in which it takes place. Obtaining information about DNA methylation at genome-wide scale has become straightforward with the use of bisulfite treatment in combination with reduced representation or whole-genome sequencing. While it is well recognized that methylation is tissue specific, a frequent limitation for many studies is that sampling-specific tissues may require sacrificing individuals, something which is generally undesirable and sometimes impossible. Instead, information about DNA methylation patterns in the blood is frequently used as a proxy tissue. This can obviously be problematic if methylation patterns in the blood do not reflect that in the relevant tissue. Understanding how, or if, DNA methylation in blood reflect DNA methylation patterns in other tissues is therefore of utmost importance if we are to make inferences about how observed differences in methylation or temporal changes in methylation can contribute to phenotypic variation. The aim of this review is to examine what we know about the potential for using blood samples in ecological epigenetic studies. I briefly outline some methods by which we can measure DNA methylation before I examine studies that have compared DNA methylation patterns across different tissues and, finally, examine how useful blood samples may be for ecological studies of DNA methylation. Ecological epigenetic studies are in their infancy, but it is paramount for the field to move forward to have detailed information about tissue and time dependence relationships in methylation to gain insights into if blood DNA methylation patterns can be a reliable bioindicator for changes in methylation that generate phenotypic variation in ecologically important traits.

Consistent scaling of inbreeding depression in space and time in a house sparrow metapopulation.

Inbreeding may increase the extinction risk of small populations. Yet, studies using modern genomic tools to investigate inbreeding depression in nature have been limited to single populations, and little is known about the dynamics of inbreeding depression in subdivided populations over time. Natural populations often experience different environmental conditions and differ in demographic history and genetic composition, characteristics that can affect the severity of inbreeding depression. We utilized extensive long-term data on more than 3,100 individuals from eight islands in an insular house sparrow metapopulation to examine the generality of inbreeding effects. Using genomic estimates of realized inbreeding, we discovered that inbred individuals had lower survival probabilities and produced fewer recruiting offspring than noninbred individuals. Inbreeding depression, measured as the decline in fitness-related traits per unit inbreeding, did not vary appreciably among populations or with time. As a consequence, populations with more resident inbreeding (due to their demographic history) paid a higher total fitness cost, evidenced by a larger variance in fitness explained by inbreeding within these populations. Our results are in contrast to the idea that effects of inbreeding generally depend on ecological factors and genetic differences among populations, and expand the understanding of inbreeding depression in natural subdivided populations.

Resistance to gapeworm parasite has both additive and dominant genetic components in house sparrows, with evolutionary consequences for ability to respond to parasite challenge.

Host-parasite relationships are likely to change over the coming decades in response to climate change and increased anthropogenic stressors. Understanding the genetic architecture of parasite resistance will aid prediction of species' responses to intensified parasite challenge. The gapeworm "Syngamus trachea" is prevalent in natural bird populations and causes symptomatic infections ranging from mild to severe. The parasite may affect ecological processes by curtailing bird populations and is important due to its propensity to spread to commercially farmed birds. Our large-scale data set on an insular house sparrow metapopulation in northern Norway includes information on gapeworm prevalence and infection intensity, allowing assessment of the genetics of parasite resistance in a natural system. To determine whether parasite resistance has a heritable genetic component, we performed variance component analyses using animal models. Resistance to gapeworm had substantial additive genetic and dominance variance, and genome-wide association studies to identify single nucleotide polymorphisms associated with gapeworm resistance yielded multiple loci linked to immune function. Together with genome partitioning results, this indicates that resistance to gapeworm is under polygenic control in the house sparrow, and probably in other bird species. Hence, our results provide the foundation needed to study any eco-evolutionary processes related to gapeworm infection, and show that it is necessary to use methods suitable for polygenic and nonadditive genetic effects on the phenotype.

Controlling for body size leads to inferential biases in the biological sciences.

Many traits correlate with body size. Studies that seek to uncover the ecological factors that drive evolutionary responses in traits typically examine these responses relative to associated changes in body size using multiple regression analysis. However, it is not well appreciated that in the presence of strongly correlated variables, the partial (i.e., relative) regression coefficients often change sign compared to the original coefficients. Such sign reversals are difficult to interpret in a biologically meaningful way, and could lead to erroneous evolutionary inferences if the true mechanism underlying the sign reversal differed from the proposed mechanism. Here, we use simulations to demonstrate that sign reversal occurs over a wide range of parameter values common in the biological sciences. Further, as a case-in-point, we review the literature on brain size evolution; a field that explores how ecological traits relate to the evolution of relative brain size (brain size relative to body size). We find that most studies show sign reversals and thus that the inferences of many studies in this field may be inconclusive. Finally, we propose some approaches to mitigating this issue.

Temporally replicated DNA methylation patterns in great tit using reduced representation bisulfite sequencing.

Seasonal timing of reproduction is an important fitness trait in many plants and animals but the underlying molecular mechanism for this trait is poorly known. DNA methylation is known to affect timing of reproduction in various organisms and is therefore a potential mechanism also in birds. Here we describe genome wide data aiming to detect temporal changes in methylation in relation to timing of breeding using artificial selection lines of great tits (Parus major) exposed to contrasting temperature treatments. Methylation levels of DNA extracted from erythrocytes were examined using reduced representation bisulfite sequencing (RRBS). In total, we obtained sequencing data from 63 libraries over four different time points from 16 birds with on average 20 million quality filtered reads per library. These data describe individual level temporal variation in DNA methylation throughout the breeding season under experimental temperature regimes and provides a resource for future studies investigating the role of temporal changes in DNA methylation in timing of reproduction.